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Abstract 

We develop a theory of the Cauchy problem for linear evolution systems of partial 
differential equations with the Caputo-Dzhrbashyan fractional derivative in the time vari- 
able t. The class of systems considered in the paper is a fractional extension of the class 
of systems of the first order in t satisfying the uniform strong parabolicity condition. 
We construct and investigate the Green matrix of the Cauchy problem. While similar 
results for the fractional diffusion equations were based on the H-function representa- 
tion of the Green matrix for equations with constant coefficients (not available in the 
general situation), here we use, as a basic tool, the subordination identity for a model 
homogeneous system. We also prove a uniqueness result based on the reduction to an 
operator-differential equation. 
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1 INTRODUCTION 



Fractional diffusion equations of the form 




(1.1) 



where < a < 1, is the Caputo-Dzhrbashyan fractional derivative, that is 




1 d_ 

r(l -a) di 



(t - r)- a u(T, x) dr - r a u(0, x) 



o 



A is a second order elliptic operator, are among the basic subjects in the theory of fractional 
differential equations. The initial motivation came from physics - the equations of the above 
type were first used for modeling anomalous diffusion on fractals by Nigmatullin [32J and for 
a description of Hamiltonian chaos by Zaslavsky [17] . See the survey papers [HI |30] for a 
description of the present status of this research area. 

The first mathematical works in this direction dealt either with the case of an abstract oper- 
ator A, that is with a kind of an abstract Cauchy problem [19] (see p2] for further references), 
or with the case where A = A is a Laplacian. For the latter case, a fundamental solution of the 
Cauchy problem (FSCP) is expressed via Fox's H-function (201111]; uniqueness theorems were 
proved in [20] for an equation with a general second order elliptic operator A; see also [37J. 
The first example of an initial-boundary value problem for this equation was considered in [46J . 
Later the initial-boundary value problems for fractional diffusion equations were studied in [28J, 
with an emphasis on the probabilistic aspects, and in [25]; for the probabilistic interpretations 
see also [i0l 129] and references therein. 

In [8] (see also [7]), Eidelman and the author constructed and investigated a FSCP for 
fractional diffusion equations with variable coefficients. We followed the classical parametrix 
method using an H-function representation for the parametrix kernel and the detailed informa- 
tion about asymptotic properties of the H-function available from [U [T7]. The fractional diffu- 
sion equation shares many essential properties with second order parabolic equations (though 
some properties are different, like, for example, the singularity of the FSCP at x = appearing 
for n>2). 

In the development of the theory of linear partial differential equations of parabolic type, 
the next step after the study of second order equations was to identify a class of systems 
which can be called parabolic. For such systems, there must be a well-posed Cauchy problem 
whose fundamental solution is an ordinary function smooth outside the singular point t = 0. 
Such a class of systems was first found by Petrowsky [31] in 1938. In fact, after Petrowsky's 
work it was understood that a complete theory of partial differential equations must include 
a thorough study of systems of equations. Petrowsky introduced and investigated not only 
parabolic systems but also hyperbolic ones, systems with correct Cauchy problems etc. For 
several decades, these subjects were central in the theory of partial differential equations; see, in 
particular, [T2l [35] . The class of parabolic systems was studied in the greatest detail [6| 17] l9~t 123] . 

So far only some special cases of partial fractional-differential systems have been considered 
[El [HJ [36l SI] . However it is the author's opinion that the fractional calculus is mature enough 
to initiate a general theory of systems of partial fractional-differential equations. In this paper 
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we follow the above line and find a fractional analog of the class of parabolic systems. Note 
that simple examples of such fractional parabolic systems appear as linearized two-component 
fractional reaction-diffusion systems used in the study of self-organization phenomena; see 
[TU| ITTj and references therein. 

We consider systems of the form (1.1) where u = (u\, . . . , u^) is a vector- valued function, 

A = A(x, D x ) = A (x, D x ) + Ai(x, D x ), (1.2) 

is a differential operator of even order 2b with matrix-valued coefficients, 

JV 

(A (x, D x )u)i = Y;Y. " ; ' {r)l) .r".r i = h---,N, (1.3) 

3=1 \P\=2b 
N 

(A 1 (x,D x )u) i = J2 4( x ) D * u i> i = h---,N, (1.4) 

3=1 \/3\<2b 

dZ = d£...D%, D Xj = -J=^-, |/?|=/3 1 + - ■• + /?„. 
We assume that all the coefficients a«(x) are bounded and satisfy the global Holder condition 

\4{x)-4{y)\<C\x-yV 

(below the letters C, c will denote various positive constants while 7 > will denote all the 
Holder exponents; for simplicity, we denote by | • | norms of all finite vectors and matrices). We 
also assume the uniform strong parabolicity condition: for all r\ e MJ 1 , z e C , 

Re(A (x,r])z,z} < -5\r]\ 2b \z\ 2 , 6 > 0. (1.5) 

In fact, the key ingredients in the construction of a FSCP for a problem with variable 
coefficients are precise estimates for the model problem 

(B { t a) U (t, z)) = My, A>(*, x) (1.6) 

containing only the homogeneous highest order differential operator in x, with "frozen" coeffi- 
cients depending on a parameter point y. As the first step, one has to consider the case where 
the coefficients clq(x), \(3\ = 2b, are constant. Already in this case, the study of a FSCP is far 
from trivial. The approach used in [HI [7] based on the H-function representation, does not work 
for systems. 

Instead, we use the subordination representation [3], [1] expressing the FSCP for the model 
system via the FSCP for the first-order (in t) parabolic system. At the first sight, it looks 
an easy approach to all fractional problems. However the subordination identity involves the 
integration in t over the half-axis (0, 00), while usually a FSCP for a parabolic equation or 
system is constructed only on a finite time interval. Nevertheless, for our model case of constant 
coefficients the subordination method works efficiently giving, by the way, new proofs of the 
estimates known for fractional diffusion equations. Note also that the probabilistic side of 
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subordination, not touched here, is also an important subject of fractional analysis; see [2j l2Tj 

SUES]. 

We also prove a uniqueness theorem for general systems (1.1). Again, the method of proving 
uniqueness in [2U1 E] (based on a kind of the maximum principle) is applicable only for second 
order equations. Here we use the reduction to an abstract equation from [JJJ] and the regularized 
resolvent estimate from [To] . 

The main results of this paper are collected in Section 2. Section 3 contains miscellaneous 
auxiliary results used subsequently. Proofs of the estimates for the estimates for the Green 
matrix of a homogeneous system with constant coefficients are given in Section 4 and are 
complemented in Section 5 with some considerations regarding the parametrix kernels. In 
Section 6, these results are used to substantiate the Levi method in our situation. The proof 
of the uniqueness theorem is given in Section 7. 

The author is grateful to the anonymous referee for helpful comments and suggestions. 



2 Main results 

In this section we introduce basic notions and formulate principal results. The proofs will be 
given in subsequent sections. 

2.1. The model system. Subordination. Let us consider systems of the form 

(®l a) u(t,xj) =A (D x )u(t,x) (2.1) 

where 

N 

{Ao{D x )u)i = J2J2 "V y >'j- i = l,---,N, 

3=1 \n\=2b 

G C, and for any rj G W\ z G C N , 

Re(A (r])z,z} < -8\r)\ 2b \z\ 2 , 5 > 0. (2.2) 

Under the assumption (2.2) (in fact, even under a weaker assumption of parabolicity in the 
sense of Petrowsky), the differential expression A (D) defines on the space L 2 (R n , C N ) of square 
integrable vector-functions with values in C N , an infinitesimal generator Aq of a Co-semigroup 
S^t) = e tAo (see [22]). 

By the subordination theorem (see Theorem 3.1 in |3J), the system (2.1) interpreted as 
an equation in L 2 (W a ,C N ), possesses a solution operator S a (t), such that for any element 

= u (x) from the domain D(Ao), the function u(t,x) = (S a (t)uo)(x), t > 0, x G M n , is a 
solution of the equation (2.1) satisfying the initial condition u(0,x) = uo{x). In addition, 



S a (t) = J tpt^S^ds, t>0, (2.3) 
o 

where (f t , a (s) = t- a $ a (st- a ), 



$ (O = V ^ 

f^ k\T(-ak + l- a y 
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so that $ Q can be written as the Wright function 



$«(C) = o*i 



(1 — a, —a 



< 



(2.4) 



See [H] for general information regarding the definition and properties of the Wright functions; 
see also Section 3.3 below. 

The function $ a is a probability density: 



OO 

$ Q (t) > 0, t>0; J <f> a {t) dt = 1. 



It is connected also with the Mittag-Lefner function 



«.(0 = Enn^. CeC ' 



fc=0 



(2.5) 



via the Laplace transform identity 

oo 

E a (-C) = J 9 a {t)e~^dt, (GC. 



(2.6) 



By the classical theory of parabolic equations, the semigroup Si(t) possesses the integral 
representation 

(S 1 (t) v )(x) = J Z(t,x-£)ri(Z)dS, V EL 2 (R n X N ), 

in terms of the FSCP Z(t,x) of the parabolic system — = Aq(D x )u. It follows from the 
estimates of $ a and Z (see below) that, for example, if rj G «S(R n ), then 

(S a (t) V )(x) = [ Z a (t,x-0v(0dC 



where 



oo 

Z a (t,x) = J ip t ,a{s)Z{s,x)ds, X^O 



(2.7) 



(as we have seen for the diffusion equations [SI Ej, Z a may have a singularity at x — 0). 
The kernel Z a is a FSCP for the system (2.1). 

In order to obtain an integral representation of a solution u(t, x) of the inhomogeneous 
equation 



it 



(t,x)j - A (D x )u(t,x) = f(t,x), u(0,x) = u Q (x), 



in the form 



t 

u(t,x) = j z a (t,x - £)«o(0 + f dT f Y ^( t - T ^ x -y)f( T ^y) d v 



(the definition of a classical solution will be given below for a more general situation), we need 
another kernel 

Y a (t,x) = (b^Z^ (t,x), x^O. 



As we will see, 



where 



Y a (t,x) = j %l)^ a {s)Z{s,x) ds, x 
o 



-sf 



(2.i 



(2.9) 



L(0,-a) 

Note that y? t a is the law of the inverse to a a-stable subordinator; see [29J. In the context 
of fractional calculus models of physical phenomena, the functions tp t ,a and ip t ,a appeared for 
the first time in [26] . 

Using properties of the functions Z and $ Q we get the integral identities 

j.a—1 



Z a (t, x) dx — 1, / Y a (t,x)dx 



f- 



r(a) 



Theorem 1. T/ie matrix-functions Z a (t,x), Y a (t,x) are infinitely differentiable for t > 0, 
i/O, and satisfy the following estimates. Denote 



R = t~ a \x\ 2 \ p(t, x) = (t- a \x\ 2b ) 2b ~ a 



(i) I}R>1, then 



D%Z a (t,x)\ < Ct~ a ^- exp(-ap(t,x)), a > 0; 



\D^Y a (t,x)\ < Ct~ a 



{n+\El 



+a-l 



exp(—ap(t, x)). 



(ii) If R< 1, n+\(3\ < 2b, then 



(iv) If R<l,n 



\D?Z a (t,x)\<Ct 
\D^Y a (t,x)\ < Ct 



,(n±l£l) 



(iii) IfR<l,n+\P\> 2b, then 



2b, then 



\D%Z a (t,x)\ < cr a \x\' n+2b ~\ p \ 
\D%Y a (t,x)\ < cr l \x\- n+2b ~w. 

\D^Z a (t,x)\ <Ct~ a , i/n = l; 



\D^Z a (t,x)\ <Cr Q [|log(r a |x| 26 )| + l], ifn>2- 



2.10) 
2.11) 

2.12) 
2.13) 

2.14) 
2.15) 

2.16) 
2.17) 



\D^Y a (t,x)\ < Cr 1 . 



(2.18) 
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(v) If R > 1, then 



(vi) If R<\, n <2b, then 



If R<1, n> 2b, then 



If R<1, n = 2b, then 



dZ a (t, x) 



dt 



<Ct 2t exp(— ap(t, x)). 



dZ a (t, x) 



dt 

dZ a (t, x) 



< Ctss 



an i 



dt 



< cr a -'\x 



a-li \-n+2b 



dZ a (t, x) 



dt 



< Cr Q - i [|log(r a |a;| 2b )| + 1] 



(2.19) 



(2.20) 



(2.21) 



(2.22) 



In all the above estimates, the constants depend only on N,n, max ajj , and the strong 
parabolicity constant 5. 



The estimates (2.10)-(2.22) agree with their counterparts for the fractional diffusion equa- 
tions [H [7], though in the latter case, for some values of n, there are more precise estimates of 
Z a and Y a . 

Note also that the fractional derivative Z a satisfies the same estimate as the derivatives 
D£Z a , \p\ = 2b. 

2.2. The general case. As stated in Introduction, we consider the system (1. 1)— (1.4) with 
bounded Holder continuous coefficients, under the uniform strong parabolicity condition (1.5). 

We call a vector-function u(t, x), < t < T, x G M n , a classical solution of the system (1.1), 
with the initial condition 

u(0,x) = uq(x), xeR n , (2.23) 

if: 

(i) u(t, x) is continuously differentiable in x up to the order 2b, for each t > 0; 

(ii) for each x G lR n , u(t,x) is continuous in t on [0, T], and its fractional integral 

t 



(Ifc a «) (t, x) 



1 



r(i-a) 



(t — t) a u(r, x) dr 



(2.24) 



is continuously differentiable in t for < t < T. 
(iii) u(t,x) satisfies the equation (1.1) and the initial condition (2.23). 
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A classical solution u(t, x) is called a uniform classical solution, if it is continuous in t 
uniformly with respect to x E R™, and the first derivative of the fractional integral (2.24) exists 
uniformly with respect to i G R". 

Our main task is to construct a Green matrix for the problem (1.1), (2.23), that is such a 
pair 

{Z^\t,x;0,Y^(t,x;0}, t E [0, T], x, £ E W 1 , 

that for any bounded function /, jointly continuous in (t, x) and locally Holder continuous in 
x uniformly with respect to t, and any bounded locally Holder continuous function -u , the 
function 

t 

u(t,x) = I Z^\t,x-t)u Q (t)dt + J d\J Y^(t-X,x;y)f(X,y)dy. (2.25) 

is a classical solution of the problem (1.1), (2. 23). 

Denote by Za\t, x — £; y) and Ya°\t, x — £; y) the kernels defined just as Z a (t, x — £) and 
Y a (t,x — £), but for the system (1.6) with the coefficients , = 2b, "frozen" at a point 
y E K™, and other coefficients set equal to zero; in (1.6), y appears as a parameter. 

Theorem 2. (a) There exists a Green matrix ^Za\t, x; £), Y&\t, x; £)j of the form 

Z<?> (t, x- = Z<°> (t, x - £ + Vz(f , x; 0, 

r« (t , x; o = rj 0) (t, x - e; + vv (t, x- , 

where the kernels Za\t,x;^), Ya°\t, x;£) satisfy the estimates listed in Theorem 1 with co- 
efficients independent on the parameter point £. The functions Vz, Vy satisfy the following 
estimates. 

(i) Ifn+\/3\ < 2b, then 

\D p x V z {t,x;£)\ < C7rt(l/ 3 l+To)| a; _^|-n+7-7o e — P(t,x-0 t o < 7o < 7, a > 0; (2.26) 

\D%V Y (t,x;£)\ < Cr 1+a -^\x - i\- n+1 e- ap{t ^ ) . (2.27) 

(ii) Ifn+\P\ > 2b, \(3\ < 2b, then 

\D%V z (t,x;Z)\ < Ct~ a+ ^\x - ^|-n+26-|/3|+ 7 -70 e -ap(M-0. ( 2 .28) 

|DfW(t,a;;0| < Cr^^lx-er"^ 26 " 1 ^ 1 ^" 27 ^-^'^. (2.29) 

(iii) // \/3\ = 2b, then 

\DlV z {t,x;i)\ < Cr a+ ^\x - i\- n+ ^e- ap ^ x - & ; (2.30) 

\D^V Y (t,x;0\ < Cr 1+fll \x - Z\- n+p2 e- ap(t ' x -° (2.31) 

where fj, 2 > 0. 
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(b) If the functions u (x), f(t,x) are bounded and globally Holder continuous (for f, uni- 
formly with respect to t), and f is continuous in t uniformly with respect to x G W 1 , then the 
solution (2.25) is a uniform classical solution. All its derivatives in x, up to the order 2b, are 
bounded and globally Holder continuous, uniformly with respect to t G [0, T]. 

Note that the estimates in Theorem 2 can be written in a variety of ways. For example, in 
(2.26) we may write 

t-&W = (t-*\z - \x - t\~W 
and obtain, taking < a' < a, that 

\D^V z {t,x;0\ < Ct-t7o| x _£|-"H/3|+7-7o e -^>*-0 
This kind of transformation is often used in proofs of various estimates in this paper. 
2.3. Uniqueness theorem. Here we maintain the same assumptions as in Theorem 2. 

Theorem 3. Let u(t,x), < t < T, x G W 1 , be a uniform classical solution of the problem 
(1.1), (2.23) with f(t, x) =0,uo(x) = 0. Suppose that the function u(t,x) and all its derivatives 
of orders < 2b are bounded and globally Holder continuous. Then u(t, x) equals zero identically. 

The rest of the paper is devoted to the proofs of the above results. Some of them can 
be extended easily to more general situations. For example, the coefficients of the subordinate 
operator A\ may depend on t. Some of the estimates of Theorem 1 (except the case n+ \/3\ = 2b) 
remain valid for systems generalizing in an obvious way the first order systems parabolic in the 
sense of Petrowsky. However for the whole range of the above results we need the condition 
(1.5). 



3 Some auxiliary results 

3.1. The Mittag-Leffter type functions of a matrix. Let B be a complex N x N matrix. The 
Mittag-Lemer function E a of the matrix B is defined by substituting B into the power series 
(2.5): 

B «( B ) = £r(TT^)' 0<a<1 - (31) 

fc=0 v ' 

Note that the analytic function E a of a matrix does not coincide with the matrix formed by 
values of the function E Q on matrix elements. 

For a class of matrices, we find a matrix analog of the asymptotic representation of the 
Mittag-Lemer function (see e.g. [5]). 

Proposition 1. Suppose that for any z G C^, 

Re(Bz,z) < -5\z\ 2 , 6 > 0. (3.2) 

Then 

E a {B) = - 1 B-' + H (3.3) 
1 (1 — a) 

where the matrix H is such that \H\ < C5~ 2 (the constant C does not depend on B). 
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Proof. Denote by 7(7", oS) the contour in the complex plane oriented in the direction of the 
increase of arg ( and consisting of the following parts: the rays 7± = {( 6 C : arg £ = ±u, \ (\> 
r} and the arc {( 6 C : — to < argC < u, |£| = r}. Here r > 0, ~ < uj < it. 

Let us use Hankel's integral representation 

' / e ? C"^C, Res>0 (3.4) 



r(s) 2717 

7(r,cj) 

(see e.g. [27]). In the integral (3.4), we make the change of variables £ = ?7 1//a . Since r > is 
arbitrary, we obtain the representation 



1 



r(s) 27T2G! 

7 (r,/3) 



e^V^- 1 ^, (3.5) 



for any /3 with ^ < (3 < rra. For our purposes, we will assume that ^ < (3 < min(|, ira), so 
that the contour 7(r, u) is located in the right half-plane, and cos £ < 0. 

Under the condition (3.2), the resolvent (B — A/) -1 exists for Re A > —5, and 



[B-XI)- 1 ]^ ' 



Re A + 5 

(see Lemma V.6.1 in [13]). In particular, \B^\ < S^ 1 . Substituting (3.4) into (3.1) we find 
that 




— [ e" 1/a (B - dri, 
na J 



7 (r,/3) 

if r > \B\. Note that on the rays contained in 7(7-, e' ?1/Q = exp (cos ^ • |77| 1//q ) where 
cos & < 0. 

a 

Let us fix r = 1. The function 

B^ [ e vl/a {B - rjl)- 1 drj 



7(1.0) 

is an analytic function of the matrix B on the open set {B : Re(I?,2, z) < — 1|^| 2 , G C Ar \{0}} 
coinciding with E a (B) on its intersection with the open set {B : |Z?| < 1} (see [21] regarding 
analytic functions on matrices; an analytic function on N x N matrices is in fact an analytic 
function of N 2 complex variables). Therefore 

E a {B) = -^- J e* 1 "* (B - rjl)- 1 drj (3.6) 

7(1.0) 
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for any matrix B satisfying (3.2). 
Next we use the identity 

(B - r]I)- 1 = B- 1 + r]B-\B - rjl)' 1 . 
Substituting it into (3.6) and using (3.5) we get 

7(1,0) 

which implies the required representation (3.3). ■ 

We will need also and estimate of E a ^ a (B) where the Mittag-Lefher type function E at(x is 
defined by the series 

Proposition 2. //a matrix B satisfies (3.2), then 

\E a , a (B)\ < C5- 2 (3.8) 

where C does not depend on B. 

Proof. As in the proof of Proposition 1, we obtain the representation (with the same 
notations) 

E a>a (B) = -^— J e^V^E-^/)- 1 ^, 

7(1,0) 

from which we get that 

E ^AB) = -^-B- 1 J ^"v'^dv-^-B- 1 J e^rf^{B - ^ drj. (3.9) 

7(1,0) 7(1,0) 

In the representation (3.5), we may pass to the limit, as s — > 0. As a result, the first integral 
in (3.9) equals 0. Estimating the second integral we come to (3.8). ■ 
In another result of this kind, we deal with the Mittag-Leffler type function 



EaAO = E 



Proposition 3. Under the assumptions of Proposition 1, 
where \H\ < C5~ 2 , and the constant does not depend on B. 
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The proof is similar to that of Proposition 1. 



3.2. Green matrices of some elliptic systems. We will need, as technical tools, estimates of 
derivatives D l3 G{x) of the Green matrix of the elliptic operator Aq(D) — I, and also estimates 
of differences D"(G(x;y') — G(x;y ; ')) of derivatives of the Green matrices of the operators 
A (y', D x )—I and A (y", D x )—I, with coefficients "frozen" at the points y' and y" . Specifically, 
we need to consider the case where n+ \/3\ = 2b. The presence of the term —I is essential - the 
matrix A Q (rj) — I (rj G M. n ) has, under the assumption (2.2), all the eigenvalues with nonzero 
real parts. For this class of elliptic systems with constant coefficients, Eidelman [6] found an 
integral representation of Green matrices, and the estimate 



\D p G(x)\ < c(l + log-^A 



\x\ < 1, (3.10) 



where it is assumed that n+|/3| = 2b (of course, other cases were considered in [6] too). It also 
follows from the constructions in [6] that G £ C 00 (R™ \ {0}). 

However the estimates for the differences are not given in [6] , and for completeness we have 
now to apply the method from [6] to this situation. 

Let N(y,t,x) be a FSCP of the parabolic system 

^^ = (A (D X )-I)u(t,x). 

Then [6] 

oo 

G{x-y) = [ N(y,t,x)dt, 



o 



so that 

£>£[G(x; y') - G(x; y")] = I D^Ntf, t, x) - N(y", t, x)} dt. 



For the FSCP N(y, t, x) we have the representation 

N(y,t,x) = (2iry n [ e^e^tf-Qt ^ = ( 27r )-«r^ e - t / e ^ 1/2b -< e My,0 d ^ 



Therefore, if n + \ f3\ = 2b, then 

Df [N(y', t, x) - N(y", t, x)} = (2n)- n r 1 j f e u ~ 1/2bx < [e Aofe '^ - e A °^'^ d£. (3.11) 

Following [6j we use the identity 



e A ( y ',ot _ e Mv",0t = I P My"m 



J e A ^"^[A (y',0 - A (y",0]e My '^dr, 
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which implies the inequality 



e Mv',t) _ e Mv",0 



< C \y'-y'Ve-W 2b . 



We can consider inequalities of this kind containing, instead of £ G MJ 1 , a point £ + ir/, 
i] G R™. By a lemma from [6] (Chapter 1), the inequality (1.5) implies the inequality 



BB(Ao(y,Z + iri)z,z) < {-5^ + M™) \ff 
(5i,/ii > 0), thus the inequality 



\ e A (y,t+i v )l < Ce -S 2 \Z\ 2b +» 2 \v\ 2b 



(62, H2 > 0). Repeating the above reasoning we get the inequality 



e A (y' £+irj) _ e A (y",Z+iv) 



<C\y'-y'Ve 



(3.12) 



{5 3 ,n 3 > 0). 

Next we use another lemma from [6] regarding the Fourier transform of an entire function 
of several complex variables satisfying an exponential estimate like (3.12). From this lemma, 
(3.12), and (3.11) we obtain the estimate 



\D^[N(y', t, x) - N(y", t, x)} I < CrVV - y'Vt 



1 26 
-Ct 26-1 Ul 26-1 



so that 



00 

\D^[G(x,y')-G(x,y")}\ <C\y'-y"[< J rV* exp ^-cT^\x\^j dt 



00 

= C\y'-y"[' J s _1 exp (-s|:r| 26 - ds = C|?/ - y"Wi{x) + I 2 {x)\ 



where 



1 1 
Ii(x) = J s~ x exp (^-s\x\ 2b - s~26^i j ds < J s" 1 exp ^-s^^rj (is < C, 



00 



,/ 2 (.r) / s _1 exp [s\x\ 2b -s-^ ) ds< I s- l e- s W 2b ds= \x\ 2b / log s • ds 
1 



00 CXJ 



This results in the required estimate 



|Df [G(ar, </) - G(ar, </')] | < C|y' - y'T (l + log , |x| < 1, n + = 2b. (3.13) 
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3.3. Properties of the functions ipt,a an d ipt,a- The function (see (2.4)) involved in the 
subordination representation (2.3) is such the § a (t) > for all t > 0, 

$ Q (t) ~ Cr l/2 e- ctTh " , t -> +oo 
(C, c > 0); see [21 US]. Changing constants if necessary, we can write the estimate 

< ipt, a {s) < ct^e-^ 1 '^ , s>0. (3.14) 
Note that (pt, a (s) — > 0, as t — > 0, for each s > 0. It is also important that 



CO 

y $ a (s)ds = 1. (3.15) 



whence 



In order to study the kernel y Q (see (2.8)), we need the function 
Making a change of variables in (2.6) we can write 

CO 

E a (-() = J ifM^ a da, 
o 

CO 

E a (-(t a ) = J n>a (a)e-^da. (3.16) 

Q 

The identity (3.16) is written in part as a warning of a complicated nature of the subordina- 
tion identities. As t — > +0, the left-hand side of (3.16) tends to 1, while in the right-hand side, 
<fit,a(cr) — > for each a > 0. Thus, it is impossible to interchange the integration and taking 
the limit. It will be reasonable here to use, instead of the Caputo-Dzhrbashyan derivative, the 
Riemann-Liouville derivative 

t 



coinciding with V> 1 ~ a u wherever u(0) = 0. In particular, if s > 0, then 

iPtM = D k°mM- ( 3 - 17 ) 

By (2.4), we have the Wright function representation 





-st~ a 


.(1 — a, —a) 
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Using the contour integral representation of the Wright function (see the equality (12.41) in 
[I]) we find that 

7+ioo 



2ni J T(l — a + aX) 

7— ioo 



s-H aA ~ a dX, 



(3.18) 



7>0, 7 ^^,^ = 0,1,2,.... 



It is known (see, for example, |18j) that D 0+ a transforms t 01 into 
Now we get from (3.17) and (3.18) that 



T(l — a + aX) 
T(aX) 



qA-1 



7+ioo 



r(A) 



T(aX) 



S -X f a\-1 dX 



7—400 



(the legitimacy of this transformation follows from the asymptotics of Gamma function). Using 
again the equivalence of the series and integral representations of the Wright functions ((12.41) 
in [1]) we prove the representaton (2.9): 



(0,-a) 



so that 



st- a ) k 



(3.19) 



(3.20) 



k\V(-ak) 

Using the asymptotics of the Wright function found in [45] (see also [4], Theorem 25) we 
get the estimate 

WtMl < Cr 1 exp {-c(sr a )T^} , (3.21) 

with c, C > 0. 

Following (2.6) and (3.16), let us write a connection of the function ipt j0t with the Mittag- 
Leffier type functions. Apply the operator -Dq+" to both sides of the equality (3.16). It is known 
(see (2.1.54) in [18J) that in the left-hand side we obtain t a ~ 1 E ctjQ (—(t a ). Thus we come to the 
identity 



(0,-a) 



st~ a 



E a , a (-Ct a ) = t- a J 0*1 


Another transformation kernel of the above kind is 
Repeating the above reasoning we find that 



e~ Cs ds, C > 0. 



(3.22) 



(3.23) 



(—a, —a) 



-sr a 



(3.24) 
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where the Wright function has the form 

l(—ct, —ol j 



k=0 



k\T(—a — ak) 



Using the asymptotics from jU H5] we obtain, for s > 0, t > 0, the estimate 

< cr^expj-^sr 01 )^} , o o. 

It is known Q18J, formula (1.10.2)) that 



(3.25) 



where E a0 (() = ~, — rr- Differentiating both sides of the identity (3.16) we find, after 
k=0 T(ak) 

elementary transformations, that 









j 0*1 




— s 


.(—a, —a) 





e-t'ds, C>0. 



(3.26) 



4 Systems with constant coefficients 

4-1. Proof of Theorem 1. Let us use the identity (2.7) which implies a representation of spatial 
derivatives, 



(4.1) 



D%Z a (t, x) = J ip tia (s)D^Z(s, x) ds, x^0. 
o 

Below we use both the integral representation for Z, 

Z{s,x) = {2ir)~ n J e ix <e sA ^) df, 

R n 

and its estimates 

I a . % | „ n+\0\ ( ,26 1 ~) 

\D p x Z(t,x)\ < Ct ~ exp j — c|x| 2b^i£~26^i I 

(c > 0) valid for all x G R n , t > (see 0). 

Note that, by its construction, the function Z a (t,x) satisfies the system (2.1) in a weak 
sense (a strong solution is obtained after a convolution in spatial variables with a test function; 
see Section 2.1). It will be seen from the investigation of the spatial derivatives below that 

A (D x )Z a exists in the classical sense for t > and x ^ 0. Therefore (dJ q) Z^j (t, x) coin- 
cides almost everywhere with A (D x )Z a (t, x) (i > 0,i ^ 0) and does not require a separate 
investigation. 



(4.2) 



(4.3) 
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It follows from the inequalities (3.14) and (4.3) that Z a (t,x) is infinitely differentiable in 
i^O, and 

oo 


(for t > 0, the second exponential factor on the right guarantees the convergence of the integral 
at infinity; for x ^ 0, the first exponential factor gives the convergence at the origin). After 
the change of variables s" 26 - 1 = a, we find that 



\DgZ a (t,x)\ < Cr a J (j 2 ^^- 26 exp {-co-|x|^} exp |-ar - ^r^=} da. (4.4) 
o 

In order to obtain an estimate for R > 1, we use the asymptotics of the integral 

oo 

11(C) = y e- ct e- dt ~"t x dt, C^oo, (4.5) 



where d > 0, x > 0, A G R, found in (formula (12.80)). Namely, 



fi(C) ~ a (dx( x )^+i exp 



:i + -)p 

X 



/T 1 / 2 (4.6) 



where p = (dx^)^, a = 2 (^-) 1/2 r(±). 

Making in (4.4) the change of variables a = t 2&-177 W e get, after easy calculations, the 
inequality 

\D£Z a {t,x)\ < Ct- a ^ j V ^ n+ ^- 2b exp (-c^) exp {-cR^-q) dr, 



^ ^ 26 1 

where the integral has the form of (4.5) with ( = i? 26 - 1 , d = c, x = , A = (n + 

1 — a 2b 

\/3\) — 2b. Using (4.6) and ignoring powers of R (changing, if necessary, the constant in the 
exponential factor) we obtain the inequality (2.10). 

Suppose that R < 1, n + \f3\ > 2b. We use again the inequality (4.4), but make the change 
of variables o = t\x\ 2b - 1 and replace the exponential factor containing \x\ by 1. We find that 



00 

\D%Z a (t,x)\ < Cr a \x\- n -M + * b J r^< 



n+\(3\)-2b e -cr dr 



2b- 1 

where — (n + \/3\) - 2b > -1. This implies (2.14). 

The inequalities (2.11) and (2.15) are proved similarly, on the basis of the estimate (3.21). 
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Let R < 1, n + \(3\ < 2b. Performing a change of variables we can rewrite (4.2) in the form 

Z(s,x) = (27t)- n s-% [ e ts ~ 1/2bx <e Ao ^ d£, 



and by virtue of (4.1), 

oo 

D?Z a (t,x) = (2n)- n r a J Q^sr^s-^r 1 ds J ^ e is ~ 1/2bx <e A °® d£, 

R" 

that is, after the change s = at a , 

oo 

= m^/..^-**/^-^* (4,) 



It follows from (2.2) that 



n + \j3\ 

(see Sect. 1.4.4 in [22]). Using the fact that $ a decays rapidly at infinity, and — — — < 1, we 

come to the inequality (2.12). The proof of (2.13), based on (2.8), (3.19), and (3.21), is similar. 

Finally, consider the most complicated case where R < 1, n + \/3\ = 2b. The representation 
(4.7) takes the form 

oo 

D?Z a (t, x) = (2iry n t~ a j <$> a {a)o~ l da J ^ e ^ a/2b ^ 1/2b -< e MO % (4 . 9) 

K« 

or, after the change £ = a 1 / 26 ^, 

oo 

D£Z a (t,x) = (2Tr)- n r a [ § a (a)do / ^e u ~ a/2bx ^e aAoM d V . (4.10) 



The identity (2.6) remains valid when the matrix —A (r]) is substituted for (. Indeed, we 
may rewrite (2.6) in the form 



oo 

k=0 ' ' { 



dt 



(the convergence of the series follows from the asymptotics of $ a and Stirling's formula). Recall 
that an entire function of a matrix is defined by substituting the matrix into the power series 
expansion. Thus we set Q = —Ao(rj), and then gather the power series into the exponential 
(keeping in mind the inequality (4.8)). Therefore 

oo 

E a (Mv)) = J ® a {v)e° A ^da. (4.11) 
o 
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If we substitute (4.11) into (4.10) and change the order of integration, we obtain the repre- 
sentation 

D?Z a (t, x) = (2n)- n t- a J rj /3 e it ~ a/2bx " n E a (A (r])) dr). (4.12) 

However this change of the order of integration requires a justification. 
It follows from (2.2) and Proposition 1 that 

\E a (A ( v ))\ < c\ v \- 2 \ \ V \>1, 

so that \r)\W \E a (A (ri))\ < C\r]\- n , \r)\ > 1. Denote 



X e (t,x,P) = (2n)- n r a / ^ a (a)a' L da / fe 



it3-it- a ' ab <r- 1 ' M x-C-A W 



e A °^> d£, e>0. 



Here we make the change of variables £ = and change the order of integration (we 

have moved away from the singularity!). Thus 



X £ (t,x,P)^(2n)- n t- a / rfe 



a / 8 it- a / 2b x-ri 



oo 

drj J <S> a (a)e aAoM da. 



Let us show that 



oo oo 

' J <S> a (a)e° A ^) d a ^ rf j <5> a (a)e° A ^da = E a (Mv))rf, (4.13) 



as functions of rj, in the topology of L 2 (lR ra ). 
Indeed, using the estimate (4.8) we see that 



$ a (a)e aAoM da 



drj 



£ 

<cJ\ V \ m Je'^da 



drj 



< J \ V \- 2n (l-e- s ^ 2b ) 2 dr, = C(I 1 + I 2 ) 



where 



h= j \r,\- 2n (l-e- s ^ 2b y dr), 



r)\-' n [l-e 



dr). 



M>i 



Using the inequality 1 — e x < x, x > 0, we find that 

h<Ce 2 j \r)\~ 2n+4b dr) = Ce 2 J \ V \ 2 ^ d V ^ 0, 



M<i 



<i 
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as e — > 0. By the dominated convergence theorem, we obtain also that I 2 — > 0, and we have 
proved (4.13). 

Now, by the properties of the Fourier transform, for any fixed t > and almost all 16K™, 
\im XJt,x,f3) = (2-K)- n r a [ ^^'^EJAoirj))^. 

£->0 J 

On the other hand, X e (t,x,p) -> D%Z a (t, x), if x 7^ 0, by (4.9). Thus we have proved the 
equality (4.12) for almost all x 7^ 0. 

Denote temporarily the right-hand side of (4.12) by X(t, x, j3). Next we prove that X(t, x, j3) 
is continuous in 1 ^ 0. This will establish the equality (4.12) for all x 7^ 0; simultaneously we 
will get the required estimate. 

By Proposition 1, X(t, x, (3) = X\ + X 2 + X 3 where 

X 1 = (2n)- n r a J V f> e«- a " a '*-*E a {A (v))dv, 
\v\<i 

f—a f 

1 P it~ '/ 2b x-r)\ 



X 2 = -— — / rfe 1 * ' ^[Mv)] d-n, 

(27r) n r(l - a) ' " ' J 



|r?| >1 



X 3 = (2ny n r a J ^e u ~ a/2bx r, H(r]) drj, 



\n\>i 

\H(r])\ < C\r]\~ 4b . Since \(3\ — Ab = —n — 2b, X 3 is continuous in x, and |Jf 3 | < Ct~ a . We see 
also that X\ is continuous in x, and \X\\ < Ct~ a . 
Let us write X 2 = X21 + X 22 + X 2 % where 

\v\>i 



(27r) n r(l - a) 
t~ a 

(27r) n r(l - a 

' \v\<i 

We have 



X ^ = -7^^ V / V^- a,2b ^[A (rj) - I]- 1 dq. 



\ [Mv) - ^\ - T+hr (4 ' 14) 

(see Lemma V.6.1 in [13]). By the resolvent identity 

[^(^)r 1 - [Mv) - 1]- 1 = [MvT^Mv) - 

we get the estimate 

20 



It follows from (4.14) and (4.15) that X21 and X23 are continuous in x, \X2i\ < Ct a , [X^l < 
Ct~ a . 

As for X22, we note that 



t 



y=t -a/2b x 



where G is the Green matrix of the elliptic operator Aq(D) — I. This means that X22 is 
continuous in x 7^ 0, thus X has the same property and coincides with D^Z a for all x 7^ 0. 
Now the required inequality (2.17) is a consequence of (3.10). 

For n = 1, we will refine this estimate. In this case, A (i]) = aorf b where 

Re(a z, z) < -5\z\ 2 for all z e C N . 

We have 

00 

D?Z a (t,x) = ^r a J ^e lt - a/2h ^E a {a r, 2b ) drj 
—00 

where 1 + (3 — 2b. In particular, the natural number (3 is odd, so that 

00 

D^Z a {t,x) = h~ a J if sm{r a/2b xr))E a {a V 2b ) drj. 


By Proposition 1, 

1 (1 — a) 

The contribution of the remainder term in the estimate of D^.Z a is clearly 0{t~ a ). Therefore 
we have to consider the function 



/ 



F (v)= f(v)sin(yr])dr], 0<y<l, (4.16) 



where cp is continuous on [0, 00), (p(rj) ~ rj 1 , rj — > 00. 

It was shown in [42J that the integral in (4.16) exists as an improper one, and 

F(y) ~ const •y" 1 ip(y' 1 ), y^O 

(asymptotics of this kind is proved in (42] for much more general situations), that is F is 
bounded near the origin. This implies the inequality (2.16). 

For the function Y a with n + \/3\ — 2b, we use, in a similar way, the identity (3.22) with a 
matrix argument £, which results in the representation 

D%Y a (t,x) = (2n)- n r 1 J ee U ~ a/2bx< E a , a (A (0)dt 
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Using Proposition 2 we obtain the desired estimate (2.18). 

d 

The estimates for the first time derivative —Z a (t, x) are obtained just as those for Z a itself. 

We use the representations (3.23), (3.24), and (3.26), as well as the estimate (3.25) and the 
matrix asymptotics given by Proposition 3. With this input, the proofs of (2.19)-(2.22) are 
similar to the ones given above. ■ 

4-2. Unified estimates. In Theorem 1, the estimates are given separately for large and small 
values of R. In order to justify the iteration procedures of the Levi method, we need unified 
estimates valid for all values of the variables. 

Proposition 4. Ifn+ \/3\ < 2b, then 

\D?Z a (t,x)\ < Cr a ^e- Cp(t ' x \ c> 0; (4.17) 

\D%Y a (t,x)\ < cr 1+a - a ^e- cp(t ' x \ (4.18) 



Ifn+ \/3\ > 2b, then 

\D%Z a (t,x)\ <Cr a \x\- n+2b -^e- cp ^ x) . (4.19) 

Ifn+ \/3\ = 2b, then 

\D%Z a (t,x)\ < Cr a [\\og(t- a \x\ 2b )\ + 1] e - cp ^ x \ (4.20) 
Ifn+ \/3\ > 2b, then 

\D%Y a (t,x)\ < Cr l \x\- n+2h -^e- cp{t > x) . (4.21) 
The constants can depend only on the parameters listed in the formulation of Theorem 1. 

Proof. The estimate (4.17) coincides with (2.10), if R > 1, being obviously equivalent to 
(2.12), if R < 1, n+ \/3\ < 2b. 

Let us consider the case where n + > 2b. It is clear that (4.19) is equivalent to (2.14), if 
R < 1. If R > 1, we rewrite the right-hand side of (2.10) as follows. Let a = d + c" (d, c" > 0). 
Then 

-c(t- a ' 2b \x\)'^=^ 



r a —e 

2b 



t~ a \x\~ n+2b ~ 



a/2b I \\ n +\P\-2b -c'{t- a / 2b \x\)-2b- a 



{t- a ' 2b \x\) 



2b 

e -c"(t- a / 2b \x\)2b- a 



< Ct^ a \x\^ n+2b ~^^e~ cp ^ t,x ^ 



where c = c", and we have proved (4.19). 

The proofs of (4.18), (4.20), and (4.21) are similar. 
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5 Parametrix 



The parametrix kernels Z a (t,x;y) and Ya (t, x; y) defined in Section 2.2 satisfy all the esti- 
mates of Theorem 1 and Proposition 4, with all the constants independent of y. 



,x;y"). These estimates are identical to those for Z a and Ya themselves, with an 
additional factor \y' — y"| 7 . The proofs are the same as in Theorem 1 and Proposition 4, with 
the following additional tools: the difference estimates for classical parabolic systems [6]; the 
difference estimate (3.13) for the Green matrices of elliptic systems; the estimate for 

e A (y',v) _ e Ao(v",v) 

given in Chapter 1 of [6]. We omit further details since they just repeat the above material. 
As in Section 2.1, we have the integral identities 



We need also estimates of the differences Z a °\t, x;y') — Z a °\t, x;y") and y^ u; (t,x;y' 



(0) 



,(0) 



Z a °\t,x;y)dx=l, / YP(t,x;y)dx 



ia— 1 



r( a y 



(5.1) 



It follows from the difference estimates and the first identity in (5.1) that 



< cr 



(5.2) 



6 The Levi method. Proof of Theorem 2 

Given the estimates of Theorem 1 and Proposition 4, the proof of Theorem 2 is carried out just 
as its counterpart for fractional diffusion equations jS] E] • The integral inequalities needed for 
the proof are given in sufficient generality in [7J. Therefore we drop the detailed calculations 
and give only the scheme and the main estimates. 

We look for the functions Z a (t, x; £), Y„ (t, x; £) appearing in Theorem 2 assuming the 
following integral representations: 

t 

Z a 1 \t,x;0 = Z a 0) (t,x-t;0 + j d\ j Y^(t - \,x - y;y)Q(X,y;£) dy; (6.1) 

K n 
t 

YW(t,x;0 = Y^(t,x-C;0 + f dxj Y^(t - \,x - y;y)$(\,y;0 dy. (6.2) 

R™ 

For the functions Q, $ we assume the integral equations 

t 

Q(t, x- = M{t, x; + j dxj K{t - A, x; y)Q(X, y; dy, (6.3) 

R n 
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t 

x- = x;0+ [d\ I K(t - A, x; y)$(A, y; (6.4) 



o 



where 

M(f, x; = [A(x, D x ) - A (t, D x )]Z®(t, x 
K(t, x; = L4(x, D x ) - MZ, D x )]Y^(t, x - £; £) 
Using the estimates from Proposition 4 we find that 

\M{t,x;0\ < Cr Q |x-er n+7 e- cp(t ' x -°, (6.5) 
|if(t,x;OI < Cr^'^x - £|-^ e -<*>(M-0 ; ( 6 . 6 ) 

c > 0, < 77 < 7. 

The increments of M and if are estimated as follows. Let A x M(t, x; £) = M(t,x;^) — 
M(£, A x K(t,x; £) = K(t,x;£) — K(t,x';^). Denote by x" one of the points x,x', for 

which \x" — £| = min{|x — £|, \x' — £|}. Then 

\A x M{t,x;0\ < Cr a \x - x'\~<~ £ \x" - £\- n+£ exp{-ap(t,x" - 0}, (6.7) 
\A x K{t,x;0\ < Cr^x - x'\~<- £ \x" - £\- n+£ exp{-ap(t,x" - 0}, (6.8) 

£ > 0. 

Using the estimates (6.5)-(6.8) we prove the convergence of iterations and obtain estimates 
for the solutions of the integral equations (6.3) and (6.4). We use, as a tool, Lemma 1.14 from 
[7J (where the conditions n = 2 and vq < 1 are in fact unnecessary). The resulting estimates 
are as follows: 

\Q(t,x;0\ <Cr a \x-i\- n ^e W {-ap{t,x-i)}, (6.9) 
|$(£,x;OI < Cr^x - £\~ n+ "< exp{-ap(t,x - (6.10) 
\A x Q(t,x;0\ < Cr a \x - x'l^lx" - ^ n+£ exp{-ap{t,x" - (6.11) 
\A x <£>(t,x;£)\ < Ct~ x \x - x'\^ e \x" - C| _n+£ exp{-ap(t,x" - £)}, (6.12) 

Now the representation of the Green matrix stated in Theorem 2 follows from its construc- 
tion (6.1)-(6.2) while the estimates (2.26)-(2.31) are obtained from (6.9)-(6.12) and Lemmas 
1.12, 1.13 from [Jj. The above estimates, together with the inequality (5.2), make it possible 
also to repeat, without significant changes, the whole reasoning from [8] or [7] regarding the 
heat potential and the initial condition. 

It also follows from (6.11), (6.12), and the difference estimates of Section 5, that the dif- 
ferences Dl [V z {t, x'; £) - V z (t, x"; £)], D% [V Y (t, x'; f) - V Y (t, x"; £)], \P\ < 26, satisfy the esti- 
mates similar to those for D@Vz, D^.V Y , with the additional factor \x' — x"| 7_e . 

Let us find out when the solution u(t, x) of the form (2.25) is a uniform classical solution 
with the global Holder properties. We consider the more complicated second term 



t 



w(t,x)= / dX / YP{t- \,x;y)f(\,y)dy = w 1 (t,x)+w 2 (t,x) 
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where 

t 

w 1 (t,x) = j dxj Y^(t-X,x-y;y)f(X,y) dy, 



o I 
t 



w 2 (t,x)= dX V Y (t- X,x;y)f(X,y)dy 



o 



(the first term in (2.25) can be considered similarly). The remainder kernel Vy is less singular 
than Ya, so that the uniform convergence of derivatives of u>2 is verified in a straightforward 
way. The same can be said about lower order derivatives of w\. 

For the leading derivatives, we have, just as in [8 J or j7], the expressions 

t 

D%w 1 (t,x) = J d\J D^°\t- X,x-y;y)[f(X,y)- f(X,x)]dy 

R n 

t 

+ J f(X,x)dX j D^°\t-X,x-y;y)dy, \p\ = 2b, (6.13) 



V) (t,x) = f(t,x) + I dX I dZ " ){t £ X y]y) [f(X,y)-f(X,x)]dy 



dt 

R n 

r(0) 



+ JnX,x)dXj 8Z ' (6.14) 

M n 

The global Holder property of the derivatives (6.13) follows from the difference estimates of 

and our assumptions regarding the function /. 
The representation (6.14) is obtained as follows (see [8] or [7j). First of all, if v(t,x) = 

(Il+ a w x ) (t,x), then B i i a) w 1 = and 

t 

v(t,x) = J dxj zW{t-X,x-y;y)f(X,y) dy. 

M« 

For a small positive number h, set 

t-h 

v h (t,x) = j dxj Z^(t-X,x-y;y)f(X,y)dy. 

R n 

dv h 



Then — — = v, + v, where 
dt 



h _ „,(1) _ 1 _„,(2) wn 

v£\t,x) = I Z^(h,x-y;y)f(t-h,y)dy, 
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wji / dA / — f(X,y)dy. 

R™ 

We have 

M = / [^(^x-^^-Z^^x-^x)] f(t-h,y) dy 

R n 

+ / 4°)(/i,x-j/;x)[/(t-/i ) y)-/(t-/i J x)]dy + /(t-/i,x). 

R™ 

It follows from the estimates of the parametrix kernel and its differences, and from the 
global Holder property of / that both integrals in the last formula tend to zero, as h — > 0, 
uniformly with respect to x G M. n , t G [0,T]. 

Similarly, we prove the convergence of v^\t, x) to the sum of the two integrals appearing in 
(6.14), uniformly with respect to x, t. This proves the uniform property of our solution. ■ 



7 Proof of Theorem 3 

First we consider the model system (2.1) with constant coefficients. A uniform classical solution 
of the system (2.1) can be interpreted as a classical solution of the operator-differential equation 

(B { t a) w) (t) = B w(t) (7.1) 

in the Banach space Cb{M. n ) N of bounded continuous vector-functions with the supremum norm. 
Here Bq is the closed operator on C&(lR n ) defined by Aq(D) with the domain 

{v G C b {R n ) N : A Q (D)v G C b (R n ) N } 

(Aq(D)v is understood in the sense of tempered distributions). Let q be such a natural number 
that q > (JV + |) • 2b. By Theorem 4.1 of the paper [16], under the condition (2.2) (in fact, the 
Petrowsky parabolicity condition would suffice), we have 

||(A/- B )-\I- A)- g/2 \\ <p(\X\), Re A > 0, (7.2) 

where p is a certain polynomial. 

Note that the operator (I—A)~ q ^ 2 is bounded on Cb(^ n ) N ] that follows from the integrability 
of its integral kernel jl3]. Under the assumption (7.2), the equation (7.1) has only a trivial 
solution w G Cfc(M n ) Ar with w(0) = 0. That is proved exactly as the uniqueness theorem from 
[T9] where it was assumed that limsup A~ 1//a log || (XI — -Bo) -1 !! = 0- O ne should just repeat 

A— >oo 

the whole reasoning from (19] for the function (I — A)~ 9//2 u>, instead of w, and notice that the 
operator (7 — A)~ g / 2 is injective. 

Thus, we have proved Theorem 3 for the system (2.1). Turning to the general case, we 
rewrite the system (1.1) with / = in the form 

(B^uj (t, x) - A (y, D x )u(t, x) = [A (x, D x ) - A (y, D x )]u{t, x) + A x {x, D x )u(t, x). (7.3) 
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Here y G M™ is an arbitrary fixed point. As before, we assume that u(0, x) = 0. 

By Theorem 2, we can write down an integral representation of a uniform classical solution 
of the Cauchy problem for the equation (7.3). By the above uniqueness result for model systems 
with constant coefficients, we obtain the equality 

t 

u(t,x) = Jdrj Y( \t-T,x-by)[MZ,Dt)-My,Dt)Ht,Z) d£ 

R n 

t 

+ JdrJ Y^(t-T,x-^y)MC,D ( )u(t,Odt (7.4) 

E" 

Denote 



X 



\m\<2b 



Differentiating both sides of (7.4) in x and taking into account the boundedness of the deriva- 
tives of u we come to the inequality 

t 

\D™u(t,x)\<C J v(r)dr J \D?Yi°\t-T,x-t;y)\ d£, \m\ < 2b. 

1" 

Using the estimates (4.18) and (4.21) we find that 

t 

\D™u{t,x)\ <C J v{r){t - r )-i+(26-H)f dTi | m | < 2b. (7.5) 
o 

For \m\ = 2b, the derivatives of a heat potential are regularized by subtraction, that is 

t 

D™u(t,x) = Jdrj D™YM(t-T,x-H;y){[Ao(t;,Ds)-A ( y ,Ds)}u(t,0 

R n 

-[A (x, D x ) - A {y, D x )]u(t, x)} d£ 

t 

Jdrj D™Y( \t-T,x-^y)lA 1 (Z,DtMt,O-A 1 (x,D x )u(t,x)]d{i. 



t 

o 



Until now, y G M. n was an arbitrary parameter. Set y = x. Then one of the terms in the first 
integral disappears. For the remaining terms, we use the inequalities 

| [Ao(Z, d*) - Ms, D t)Mt, 01 <c\x- t\Mt), 

\A^, D ( )u(t,0 - A l {x,D x )u{t,x)\ < C\x - £| sup \D^ x u(t,x)\ < C\x - £\v(t) 



x€R n 



(since degAi + 1 < 2b). 
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It follows from (4.21) that 

t 

\D™u(t,x)\<C J v{r){t - t)- 1 dr 



o 

X 



J (\x - er + \x- ei)k - er n ex P {- c [(* - r)-> - ^ ^ } % 

t™ 

< Cy w(r)(t-r)- 1+ ^dr, 
o 

|m| = 26. Adding this inequality to (7.5) we find that, for any t G [0, T], 

t 

v(t) <C I v(r)(t - t)- 1+£ dr, e > 0. (7.6) 



By a kind of the Bellman- Gronwall inequality proved by Henry (|15j. Lemma 7.1.1), it follows 
from (7.6) that v(t) = 0, hence u(t, x) = 0. ■ 
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